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Abstract 

m 

•^ , We analyze the most salient cosmological features of axions in extensions of the Standard Model 

tH- I with a gauged anomalous extra U{1) symmetry. The model is built by imposing the constraint of 

['nT ' gauge invariance in the anomalous effective action, which is extended with Wess-Zumino countert- 

erms. These generate axion-like interactions of the axions to the gauge fields and a gauged shift 

symmetry. The scalar sector is assumed to acquire a non-perturbative potential after inflation, at 

^^ , the electroweak phase transition, which induces a mixing of the Stiickelberg field of the model with 

the scalars of the electroweak sector, and at the QCD phase transition. We discuss the possible 

mechanisms of sequential misalignments which could affect the axions of these models, and gen- 

j^ ' erated, in this case, at both transitions. We compute the contribution of these particles to dark 

^ . matter, quantifying their relic densities as a function of the Stiickelberg mass. We also show that 

models with a single anomalous U(l) in general do not account for the dark energy, due to the 

presence of mixed U{1) — SU{3) anomalies. 
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1 Introduction 

Given its important role as a possible solution of the strong CP problem llH as well as a candidate 
for the dark matter of .he urriverse, the study of axiorrs 0, 1 0, g fi B (- i for arr overview) has 
received momentum both at theoretical and experimental level along the years. The invisible axion 
owes its origin to a global U{1)pq (Peccei-Quinn, PQ) symmetry which is spontaneously broken in the 
early universe and explicitly broken to a discrete Z^ symmetry by instanton effects at the QCD phase 
transition [9|. The breaking occurs at a temperature Tpq below which the symmetry is nonlinearly 
realized. Strings and domain walls relics, which are typical of axion models and are a problem in 
ordinary PQ cosmology, can be avoided by introducing inflation to account for their dilution, or by 
embedding the model into more general constructions based on theories of Grand Unification [lOt]. 

The almost massless nature of the axion and its suppressed coupling to the fields of the Standard 
Model are consequences of the fact that this field is associated with the phase of a global anomalous 
symmetry. Both properties are related to the same scale, the axion decay constant fa ~ 10^'^ — 10^^ 
GeV. 

The implications of the PQ axion in cosmology, both in supersymmetric and in non supersymmetric 
models, have been explored to a finer level of detail. For instance, the axion plays an important role in 
determining the structure of the primordial perturbations [lll.ll2l. Il3l|. where it can act as a curvaton. 

The gauging of an anomalous symmetry has some important effects on the properties of this 
pseudoscalar, first among all the appearance of independent mass and couplings to the gauge fields. 
This scenario allows a wider region of parameter space where to look for these particles. For this 
reason, axion-like fields, which are at the center of several investigations, are unlikely to find any 
significant and fundamental formulation without an underlying anomalous U{1) gauge symmetry, as 
emphasized in previous works [ij, ll5[ , [16|] . 

So far only two complete models have bee n p ut forward for a consistent analysis of these types of 
particles, the MLSOM [ITj] and the USSM-A []J]. The first of them is at the basis of the elaborations 
that we are going to provide in this work. Here we will be focusing on the phenomenological analysis 
of a scenario which is a direct consequence of the model introduced in [17], while more details on the 
supersymmetric construction will be discussed in a separate work. 

Although the natural framework that motivates these constructions is open string theory [l9| , the 
effective actions describing these types of particles can be consistently defined at lower energy just 
by the inclusion of the relevant dimension-5 Wess-Zumino (Peccei-Quinn) interactions. These are 
necessary in order to guarantee the gauge invariance of the effective action and can be interpreted as 
counterterms. In fact, they balance the anomalous variation of the 1-loop effective action induced by 
the extra U{1) symmetry, restoring the gauge symmetry. 

The gauging of an anomalous symmetry is the essential element in the construction of these 
effective actions and can be justified within intersecting brane models. The gauging is a variant of the 
standard Peccei-Quinn construction and is characterized by a new scale M, which is the Stiickelberg 



mass. We recall that Stiickelberg extensions of the Standard Model with a non-anomalous U{1) have 



been analyzed in several recent works 20. l2ll. |22|. 



We are going to provide a physical perspective on the possible phenomenological implications of 
the anomalous case. In particular, we will try to connect the Stiickelberg fields, which are in the 
spectrum of these models, to the physical axion which may appear as an extremely weakly interacting 
particle of a certain relic density in our current universe. Our assumption, in the identification of 
the physical axion, is that the original Stiickelberg fields will mix at the electroweak phase transition 
with the Higgs sector. As a result, an almost massless state will emerge after the electroweak phase 
transition. 

One of the key mechanisms that we will try to adapt and extend from the PQ case is that of vacuum 
misalignment. This phenomenon occurs whenever a quasi Nambu-Goldstone mode - generated by the 
breaking of a certain symmetry - acquires non-perturbatively a small potential, lifting one flat direction 
from the vacuum degeneracy. For axions characterized both by an SU{2) and an SU{3) charge the 
mechanism of vacuum misalignment becomes sequential, as we are going to show. 

From a more general perspective, we will also try to characterize the possible role of these types of 
particles as quintessence axions. These appear in models where the axions remain decoupled from the 
gluonic sector and their mass is purely of electroweak origin (see for instance [23|). In this case one 
tries to exploit the Nambu-Goldstone nature of these particles. The main idea behind this proposal is 
that a phase transition around the electroweak scale can generate a small curvature in the potential, 
capable of giving a tiny mass to this particle, smaller than the Hubble parameter at current time 
(Hq). For this to be possible, as we are going to show, one has to search for solutions of the anomaly 
equations for an anomalous U{1) which has a vanishing mixed anomaly with the SU{3) color group. 
In this case the only source of mass for these axions would come from the electroweak and not from 
the QCD phase transition, and as such could be extremely small. 

In the general models that we analyze, the anomaly equations do not allow for such a solution, 
although this would not exclude the possibility of finding others, in the presence of more complicated 
gauge structures, for instance in models with several U(l)'s. We will not address this specific point 
any further, leaving it as an option for future studies. Instead, we will concentrate on the general 
features of an axion-like field coming from a single anomalous U{1) symmetry, characterized by the 
presence of mixed anomalies both with the SU{2) and SU{3) sectors. The phenomenological details of 
the model are rather intricate, and have been worked out before. For this reason we have summarized 
in the next section some of their salient features, which turn out to be necessary in order to proceed 
with a realistic estimate of the relic densities. This is the specific goal of our work. 



2 General features of models with gauged axions: the Stiickelberg 
field 

In this section we briefly review the main features of the class of models that we address, discussing 
specifically the Stiickelberg field b which accompanies their anomalous U{1)b symmetry. It has been 
included in order to clarify the origin of the anomalous gauging and to compare the roles played by the 
PQ (a) and the Stiickelberg axions, which is relevant for the analysis that will follow. The structure of 
the entire Lagrangian is discussed in 17| and has been briefly summarized, in part, in the appendix. 



Intersecting brane models are one of those constructions where these types of generalized axions 



appear 
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25l . |26| |. In the case in which several stacks of branes are introduced, each stack being 
the domain in which fields with the gauge symmetry U{N) live, several intersecting stacks generate, 
at their common intersections, fields with the quantum numbers of all the unitary gauge groups 
of the construction, such as U{Ni) x [/(iVa) x ... x U{Nk) = SU{Ni) x [/(I) x SU{N2) x C/(l) x 
... X SU{Nk) X U{1). In realistic models, the phases of the extra C/(l)'s are rearranged in terms of 
an anomaly-free generator, with an (anomaly free) hypercharge U(l) (or C/(l)y) times extra C^(l)'s 
which are anomalous, carrying both their own anomalies and the mixed anomalies with all the fields 
of the Standard Model. 

For instance, a simple realization of the Standard Model is obtained by taking 3 stacks of branes: 
a first stack of 3 branes, with a symmetry f/(3), a second stack of 2 branes, with a symmetry U{2) and 
an extra single brane U{1), giving a gauge structure of the form SU{3) x SU{2) x [/(I) x [/(I) x C/(l). 
Linear combinations of the generators of the three C/(l)'s allow to rewrite the entire abelian symmetry 
in the form U{1)y x U{iy x [/(I)". These rearrangements of the U{1) phases have been studied in 
the previous literature. For instance, the original basis for the U(l)'s is also called "the brane basis", 
while the reorganization of the generators in the form of "hypercharge plus reminder" goes under the 



name of "the hypercharge basis" . There are explicit assignments in the recent literature [2J, |25|, [27 1 . 
The simplest realization of the Standard Models (SM) is obtained by 2 stacks and a single brane at 
their intersections, giving a symmetry C/(3) x U{2) x C/(l). In this case, in the hypercharge basis, the 
gauge structure of the model can be rewritten in the form SU{3)c x SU{2)w x C/(l)y x ?7(1)' x U{1)" . 
We will be using also the notation U{1)b x U{1)c to refer to the two C/(l) factors (f/(l)' x U{1)") 
of the abelian gauge structure. As often emphasized in previous works, the two extra U(l)'s are in a 
"broken" phase. For instance, if we denote with B and C, the kinetic terms of these abelian anomalous 
gauge fields are given by 






^St = 7T (9^6 - MiB^Y + - {d^c - M2C^y , (1) 



which is the well-known Stiickelberg form. Mi and M2 are also called Stiickelberg masses while b 
and c are two pseudoscalars known as Stiickelberg fields (or "Stiickelberg axions"). The Stiickelberg 
symmetry of the Lagrangian ([1]) is revealed by acting with gauge transformations of the gauge fields 



B and C, under which their corresponding axions b and c vary by a local shift 

SbB^ = d^OB Sb = Mi9b 

5cC^ = d^9c 6c = M29c, 

parameterized by the local gauge parameters 6b and 6c- In the literature, the Stiickelberg symmetry 
is presented as a way to give a mass to an abelian gauge field but still preserving the gauge symmetry 
of the theory. However, a more careful look at this symmetry shows that its realization is the same 
one obtained, for instance, in an abelian Higgs model when one decouples the radial excitations of the 
Higgs fields from its phase [IGj]. Therefore, in this respect, the symmetry does not appear to contain 
much novelty. However, in the effective theory which characterizes these models, the mechanism which 
generates the mass of the anomalous [/(l)'s is unrelated to the traditional Higgs mechanism, since 
there is no Higgs potential involved. 

The massive anomalous gauge bosons acquire a mass through the presence of "^ A F " couplings in 
the effective string theory description (see for instance ^8|] ) . The starting Lagrangian of the effective 
theory involves an antisymmetric rank-2 tensor A^^j coupled to the field strength F^j^ of an anomalous 
gauge boson (here denoted by B) 

^2^ ^t^up ^^2^ ^t,u+ ^ 

where 

J^fiup = O^Ayp + OpA^y + OpApp^, r ^i, = Ofj^Bn — OuBf^ [S) 

is the kinetic term for the 2-form and g is an arbitrary constant. Beside the two kinetic terms for Ap_i, 
and Bp^, the third contribution in Eq. ([2]) is the A A F interaction. 

The Lagrangian is dualized by using a "first order" formalism, where H is treated independently 
from the antisymmetric field Apy. This is obtained by introducing a constraint with a Lagrangian 
multiplier field b{x) in order to enforce the condition H = dA from the equations of motion of b, in 
the form 

_fJIJ-!^PfJ _PP-'^ P _ _ fP-^P<^fJ B +- 

-^2 J^fiup 4„2 P" g ^ -n^i^p -Do- -I- g 



^ = -^H'^'^'Hpyp - -^F^^'^Fp, + — e^^'^P'^Ap, Fp^, (2) 



£o = --HP'^PRpyp - —F'^'' Fp, - — e^'^P'^Hp^p B, + - b{x) e^'^P'^dpH.p^. (4) 



The appearance of a scale M in this Lagrangian is of paramount importance both in the analysis 
of the relic densities of axions generated by the dualization of this action, and in determining the 
mass of the extra anomalous U(l) gauge boson, which has been analyzed in detail in previous works 



29|. It defines the energy region where the Green-Schwarz mechanism comes into play to cancel the 



anomaly in orientifold vacua of string theory [17|] . Clearly, it is part of a far more involved field theory 
Lagrangian which, in general, is not included in the field theory analysis of this mechanism, since the 
expansion stops at operators of dimension 5. We just remark, at this point, that the appearance of 
the Stiickelberg description in theories with gauge anomalies is not limited to effective field theories 



derived from strings, but it is also common to simple 2-dimensional models, such as the bosonized 
Schwinger model (see [30]). 

The last term in ^ is necessary in order to reobtain ([2]) from (j3]). If, instead, we integrate by 
parts the last term of the Lagrangian given in (j3]) and solve trivially for H we find 

Hf^'P = -eP'^P" {MB„ - d^h) . (5) 

Inserting this back into ^ we obtain the expression 



1 1 

p^^^ /? _ ± 



Ca = -^r^F^^ F^,--{MB^-d„bf (6) 



which is the Stiickelberg form for the mass terms of B. 
This rearrangement of the degrees of freedom, valid in a classical sense [31|, and the mapping of 
the possible physical phases of these two model theories, is an example of the connection between 
Lagrangians of antisymmetric tensor fields and their dual formulations, that in this specific case is an 
abelian massive Yang- Mills theory in a Stiickelberg form (see for instance the discussion in [S^]). 

The axion field generated by the dualization mechanism appears to be a Nambu-Goldstone mode, 
which could be absorbed by a unitary gauge choice in the (defining) Stiickelberg phase of the model. 
However, as discussed in [l7[, we will allow a mixing between this mode and the Higgs sector at the 
electroweak phase transition, by introducing an extra potential which respects the gauge symmetry 
and whose origin has been left, so far, unspecified. This mixing potential is here assumed to be of 
non-perturbative origin and triggered at the electroweak phase transition. It is parameterized by 
constants (Aj) which are strongly suppressed by the exponential factor (~ g""^'""*, with Sinst the 
instanton action) , determined by the value of the action on the instanton background (for electroweak 
instantons). We will come to discuss these points rather closely in the next sections. 

For this reason, at low energy, the counting of the physical degrees of freedom in the pseudoscalar 
sector of the model is performed in the combined Higgs-Stiickelberg phase, where a massive physical 
axion emerges from the combination of the phases of the Higgses and of the Stiickelberg field. In 
models with several C/(l)'s this construction is slightly more involved, but the result of the mixing of 



the complex CP odd phases leaves as a remnant, also in this case, a physical axion, denoted by x 13], 
whose mass is controlled by the size of the Higgs-axion mixing. 

The Stiickelberg Lagrangian that we have reviewed is part of the classical action 5*0 which also 
includes the remaining gauge kinetic terms of the theory at classical level, for a symmetry SU{'i) x 
SU{2) X [/(l)y. The remaining interactions can be found in the appendix. 

2.1 Charge assignments and counterterms 

We refer to the appendix for more details concerning this class of models and for our conventions, 
together with a brief outline of the structure of the counterterms in the effective Lagrangian. Here we 
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Figure 1; Anomalous contributions to the Lagrangian and WZ eountertcrms 



briefly comment on the list of the charge assignments of the single extra U(l) model, which is given 
in Table ([I]). 

Specifically, q^^qn denote the charges of the left-handed lepton doublet (L) and of the quark 
doublet (Q), while 5^^, g^, g^ are the charges of the right-handed SU{2) singlets (quarks and leptons). 
We denote with Ag-^ = qS ~ Qd ^^^ difference between the two charges of the up and down Higgses 
{Qu J qd) respectively. The trilinear anomalous gauge interactions induced by the anomalous f7(l) and 
the relative counterterms, which are all parts of the 1-loop effective action, are illustrated in Fig. [TJ 
The numerical values of the counterterms appearing on the second line of Fig. [1] are fixed by the 
conditions of gauge invariance of the Lagrangian and are summarized by the following relations 
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CBgg = -i-2q^ + qS^+qS^), 



1 



Cbww = 7;{-1l - Mq) 



(7) 



They are, respectively, the counterterms for the cancellation of the mixed anomaly U{\)bU{1)\ and 
[/(l)y[/(l)^; the counterterm for the BBB anomaly vertex or U{\)\ anomaly, and those of the 
U{1)bSU{?>)'^ and U{1)bSU{2)'^ anomalies. They are defined in the appendix. From the Yukawa 
couplings we get the following constraints on the U{\)b charges 



QQ-qd-ldn 



Qq+Qu 



dun 



^L 



l^d 



7^ 



(8) 



In Tab. ([T]) we also show the expressions of the free U{\)b charges appearing on each generation, having 



taken into account the conditions of gauge invariance of the Yukawa couphngs. Using the equations 
above, we can ehminate some of the charges in the expression of the counterterms, obtaining 

1 



CBYY = -i3qf+9q^ + 8Aq^), 

Cybb = 2 [qfiqE + 3g|) + 2Aq^{qS + q^) + {Aq^f] , 

Cbbb = (gf - q^f + 3(<zf + q^ + Aqy + 3{q^ - q^f - 2{qff - 6{q^f, 
Aq^ 



CBWW = \i-qf-3q^)- (9) 



The solutions given above are generic, in the sense that they parameterize, in principle, an infinite class 
of models whose charge assignments under U{1)b are arbitrary, with the charges on the last column 
of Tab. ([l]) taken as their free parameters. One can immediately observe that, due to the presence, in 
general, of a nonvanishing mixed anomaly of the U{1)b with both SU{2) and SU{3), the Stiickelberg 
axion of the model has interactions both with the strong and the weak sectors, which support instanton 
solutions, and therefore could acquire a mass non-perturbatively both at the electroweak and at the 
QCD phase transitions. Notice, in particular, that for a model in which Aq = 0, in which both 
doublets of the Higgs sector, H^ and Hii carry the same charge under U{1)b, then the axion mass will 
not acquire any instanton correction at the QCD phase transition. At the same time, however, it is 
easy to show that in this case the potential responsible for Higgs-axion mixing disappears. Therefore 
the axion remains a Nambu-Goldstone mode which is completely absorbed at the electroweak phase 
transition. In this case, obviously, there is no mechanism of vacuum misalignment for the axion field 
(b). This will contribute to the mass of the two neutral gauge bosons Z and Z' , just like all the neutral 
components of the two Higgses of the model. 

The solution of the same equations with a vanishing electroweak interactions of the Stiickelberg 
appears instead possible by choosing q^ = —^qq. In the presence of both a weak {Cbww) and a 
strong (CBgg) counterterm, we will assume that the massless Stiickelberg field b will mix with the 
scalar CP-odd sector and a physical axion (x) will emerge from this mixing with a tiny mass (m^) 
generated by electroweak instantons. The corresponding potential will be rather shallow and for this 
reason this new degrees of freedom will be essentially misaligned but frozen. Its contribution to the 
relic density will be indeed negligible and for this reason at this stage x is extremely light, and massless 
for all practical purposes. However, due to the presence of a coupling of this field with the strong 
sector, its mass will be significantly modified at the QCD phase transition, as in the Peccei-Quinn 
case, with a value which will depend on the size of the Stiickelberg mass M. 
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Table 1 : Charges of the fermion and of the scalar fields 



3 The electroweak potential for massless fields 



As in previous works [15*], in the construction of the effective action we follow a bottom-up approach 
with general charge assignments parameterized just by the set of free charges of \J{\)b- These are 
shown in Fig. [H together with the fundamental gauge structure of the Standard Model. The scalar 
sector of the anomalous abelian models that we are interested in is characterized by a rather standard 
electroweak potential involving, in the simplest formulation, two Higgs doublets Vpq{Hu, H^) plus 
one extra contribution, denoted as Vf>0{Hu,Hd,b) or V' , [l7| which mixes the Higgs sector with the 
Stiickelberg axion b, needed for the restoration of the gauge invariance of the effective Lagrangian 



V = Vpq{H^, Ha) + V>m(F„, Hd, b). 



(10) 



The appearance of the physical axion in the spectrum of the model takes place after that the phase- 
dependent terms, here assumed to be of non-perturbative origin and generated at the electroweak 
phase transition, find their way in the dynamics of the model and induce a curvature on the scalar 
potential. The mixing induced in the CP-odd sector determines the presence of a linear combination 
of the Stiickelberg field b and of the Goldstones of the CP-odd sector, called x-, which is characterized 
by an almost flat direction. To better illustrate this point, we begin our analysis by turning to the 
ordinary potential of 2 Higgs doublets. 



VPQ 



tJ-uHlHu 



f^lHlH, 



XuuiHiH^f + \dd{HlHdf - 2\UHiHu){H\Hd) + 2X'^a\H^T2Hd\'' 

(11) 



to which we add a second term 

UHlHa){HiHae-'<^B^'i--'i^)^) + h.c., (12) 

These terms are ahowed by the symmetry of the model and are parameterized by one dimensionful (Aq) 
and three dimensionless constants (Ai,A2,A3). They are assumed to be generated at the electroweak 
phase transition non-perturbatively, and as such their values are related to an exponential factor 
containing as a suppression the instanton action. In the equations below we will rescale Aq by the 
electroweak scale v = \ v'^ + "u^ (Aq = Aof ) so to obtain a homogeneous expression of the mass of x as 
a function of the relevant scales of the model which are, beside the electroweak vev f , the Stiickelberg 
mass M and the anomalous gauge coupling of the U{1)b-, gs- 

The physical axion x emerges as a linear combination of the phases of the various terms, which 
are either due to the components of the Higgs sector or to the Stiickelberg field b. To illustrate the 
appearance of a physical direction in the phase of the extra potential, we focus our attention just on 
the CP-odd sector of the total potential, which is the only one that is relevant for our discussion. The 
expansion of this potential around the electroweak vacuum is given by the parameterization 

This potential is characterized by two null eigenvalues corresponding to two neutral Goldstone modes 
{Gq,Gq) and an eigenvalue corresponding to a massive state with an axion component (x). In the 
(lm.H^,lm.H^,b) CP-odd basis we get the following normalized eigenstates 

Gl = I {vd,Vu,^) 



^l- . ^ i 9B{<l,-<lu)vavl9B{,,-qu)vlvur...L.,..2^ 



- / , = / — 1 / — ' V ^ivi y ^^u^ ^d 



'K + Vd \n + < 



X 



V2Mvu, -V2Mvd,gB{qd - qu)vdVu) (14) 



/g%iqd-qu?vlvj + 2M^vj + vl) 
and we indicate with O^ the orthogonal matrix which allows to rotate them on the physical basis 



Gq 



^0 

X 




Gl\=0^\ Imi^o I , (15) 
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which is given by 



O 



X 



f Ha i^ \ 

gB{qd-qu)vd-vl 9B(qd-Qu)vlvv. y/2Mv 

v^ g%{qci~qu)^vlvl+2M'^v2 v^ g%(qd~qufvlvl+2M^v^ yj g%(qd~qu?vlvl+2M'^v'^ 

. 'j2Mvu V^Mvg 9B{qd-qu)vdVu 

\ V9Bild-qu)^v'ivl+2M2v2 ^/gl{qd~qu)^vlvl+2M^v^ ^ 9\{qd~qu?vlvl+2M-^v^ ) 



(16) 



where f = a /w^ + "V^- 

X inherits WZ interaction since h can be related to the physical axion x ^-iid to the Goldstone modes 

via this matrix 

h = O^Gl + O^Gl + Olx, (17) 

or, conversely, 

X = 0^^1mHd + 0^2^mHu + 0^^b. (18) 

Notice that the rotation of b into the physical axion x involves a factor O33 which is of order v/M. This 
carries as a consequence that x inherits from b an interaction with the gauge fields which is suppressed 
by a scale M'^/v. This scale is the product of two contributions: a 1/M suppression coming from the 
original Wess-Zumino counterterm of the Lagrangian {b/MFF) and a factor v/M obtained by the 
projection of b into x due to O^. 

More details on the structure of the various operators appearing in this model have been included 
in an appendix in order to make our treatment self-contained. We have included also a brief discussion 
of the construction of (7x77' '^hich is the factor in front of one of the most important counterterms 
needed in our numerical analysis and which controls the decay of the axion into photons. We briefly 
comment on its structure. 

The final coupling appears as a coefficient in the interaction of the physical axion with two photons 

5x77^ 7 7 l-'-^J 

and is given by 

g^^ = (FOil.3^0^3, + CyyO^^O^^) 0^3. (20) 

It is defined by a combination of matrix elements of the rotation matrices O^ and O^, together with 
some counterterms F and Cyy- O is the matrix that rotates the neutral gauge bosons from the 
interaction to the mass eigenstates after electroweak symmetry breaking and has elements which are 
0(1), being expressed in terms of ratios of coupling constants. They correspond to mixing angles. 
The coefficients F and Cyy are the W Z counterterms for cancelling the anomalies emerging from the 
SU{2)U{1)\ and U{1)bU{1)\ sectors and can be found in the appendix. They are both suppressed 
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by 1/M, while the matrix element O33, as we have mentioned, is of order v/M. Defining g^ = §2 + Qy: 
the expression of this coefficient can be given in the form 



^2 „2 



"- = IS0O- T. {-"?'■ + "h (lUf - ifL U].)') ■ (21) 

Notice that this expression is cubic in the gauge coupling constants, since factors such as 52/5 and 
qy / 9 are mixing angles while the factor I/tt^ originates from the anomaly. Therefore one obtains a 
general behaviour for g^^ of 0{g^v/M'^), with charges which are, in general, of order unity. 

3.1 Periodicity of the V potential 

The phase-dependent potential has a well-defined periodicity. To identify the corresponding phase in 
the Higgs-neutral CP-odd sector we introduce a polar parametrization of the neutral components in 
the broken electroweak phase 

i^o = -^ (^V2vu + pUx)) e^^ H^a = -^ (fiv^ + />|l(x)) e*^^, (22) 

where we have introduced the two phases -F„ and -F^^ of the two neutral Higgs fields. The potential is 
periodic with respect to the linear combination of fields 

m - ^^%^K-) - -^Pl{^) + -^Fl{-\ (23) 

2M ^J2vu V^Vd 

and using the matrix O^ to rotate on the physical basis, the phase describing the periodicity of the 
potential turns out to be proportional to the physical axion, modulo a dimensionful constant (a^) 

e{x) = ^, (24) 

where we have defined 

2vuVdM 
- - (25) 



/gUQ<i-qurvjvl + 2M^ivj+vl) 

Notice that a^, in our case, takes the role of fa of the PQ case, where the angle of misalignment is 
identified by the ratio a/ fa, with a the PQ axion. In our case ay., however, is of the order of the 
electroweak scale. This, as we are going to show, has drastic implications on the relic densities of 
axions generated at this transition. 
Notice that x (o^j equivalently, 6) is gauge invariant as one can check quite directly. In fact a U{1)b 
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infinitesimal gauge transformation with gauge parameter aB{x) gives 

i 



5Hu = 


' --^qugBOBHu 


SHd = 


-- —^qdOBasHd 


6F- = 


Vu 

qugBOB 


6F,^ = 


Vd 
prfldgBOLB 



6b = -Mas (26) 

giving 69 = 0. The gauge invariance under U{1)y can be easily proven by using the invariance of 
the Stiickelberg field b and the fact that the hypercharges of the two Higgses are equal. Finally, the 
invariance under SU{2) is obvious since the linear combination of the phases that define 6{x) are not 
touched by the transformation. From the Peccei-Quinn breaking potential we can extract the following 
periodic potential 

V =4vuVd {X2vj + Xsvl + Ao) cos (^^ + 2Xivlvj cos ^2^^ , (27) 

with a mass for the physical axion x given by 



m 



2 _ 2VuVd ^^ _2 I \ „,2 , \ „,2 , ^ \ ,, „, \ ^ \„,2 

X 2 

"x 



{Xov^ + X2VJ + A3< + 4.XiVuVd) ^ Xv\ (28) 



Notice that, according to our assumption about the origin of the extra potential, this is driven 
by the combined product of non-perturbative effects, due to the exponentially small parameters 
(Ao, Ai, A2, A3), with the electroweak vevs of the two Higgses. Notice also the irrelevance of the 
Stiickelberg scale M in determining the value of a^ ~ 0{v) and of m^ near the transition region, 
due to the large suppression factor A in Eq. (f28]l . One point that needs to be stressed is the fact that 
at the electroweak epoch the angle of misalignment generated by the extra potential is parameterized 
by x/<^x while the interaction of the physical axion with the gauge fields is suppressed by M'^/v. 
This feature is obviously unusual, since in the PQ case both scales are a single scale, the axion decay 
constant fa- 

We will consider in the next sections two possible scenarios, the first is the low gravity scenario, 
where M^/w is in the TeV region or above, but essentially disconnected from the typical scale appearing 
in typical PQ axion models {fa ~ 10^^ GeV). In the second scenario we will allow a very large value 
for M^/f , of the same order of fa- In this second case we will re-obtain the PQ axion model, with 
relic densities for x which are comparable with those typical of a PQ axion. The appearance of 
a misalignment angle of the form xl^x-i respect to the PQ case {a/ fa), is going to have drastic 
consequences on the relic densities of this particle generated at the electroweak scale, densities which 
will be found to be negligible. 
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At the QCD phase transition a new - much more sizeable - misahgnment occurs and the axion 
mass gets enhanced by the QCD instantons respect to m^ given in (I28p . which in this case is typically 
of 0{Aq(-,j^v/M'^). Notice that at the QCD phase transition M'^/v takes the same role of fa in PQ 
(in the PQ case m^ ~ ^Qco/fa))- 

There are also some crucial points of difference between a gauged axion and the PQ case that require 
some comment, since they are not so obvious. Notice, in fact, that in the PQ case, if vpQ is larger 
than the scale of inflation, then the value of the 6 field can be considered essentially homogeneous. 
We have already mentioned that in the PQ case is a physical field at every physical scale, since it 
is a Nambu-Goldstone mode of a global U{1) symmetry and as such cannot be gauged away. In that 
case the role of the mechanism of vacuum misalignment at the QCD transition is just to provide a 
mass for this Goldstone mode. 

In our case, instead, b has no potential and is charged under an anomalous gauge symmetry. As 
such it appears as a longitudinal component of the anomalous gauge boson B, above the electroweak 
scale. This also implies that there is no effect on b due to infiation, being the Stiickelberg not a 
physical field at the scale of infiation. Thus, the reappearance in the CP-odd sector of a component 
of 6 as a physical axion, x, at the electroweak phase transition, implies that this physical component 
is not a homogeneous field at the electroweak time. Similar types of inhomogeneities are found also 
in the PQ case, in models characterized by late infiation. In fact, in that case the homogeneity of the 
axion field beyond the QCD horizon is not guarantee either. 

For this reason, it is conceivable that x is homogeneous within the electroweak horizon for the 
same argument, and one can neglect fluctuations of x that enter the horizon at later times. Anyhow, 
even if these fluctuations were included, they are likely to play a minor role respect to other, more 
significant effects, such as those determined by the size of the Stiickelberg mass M (M'^/v), which 
has a dominant impact on the value of the relic densities for these types of axions. For this reason 
we will be leaving aside possible further corrections due to a non-homogeneity of the b field beyond 
the electroweak horizon, knowing that variants of this approach could be worked out following the 
discussion given in (33i |. 



4 Decays of axion-like particles 

The physical axion acquires a non- vanishing coupling with the massive fermions that is proportional to 
the rotation matrix O^ and to the mass of the fermion. This coupling increases the number of its decay 
modes and, in particular, induces new channels in its decay rate into gauge bosons, mediated both by 
fermion loops and by direct Wess-Zumino interactions. In this section we perform a complete study of 
the decay rate under the assumption that the mass of x is in the meV region and below. In particular, 
in the case of a very light axi-Higgs, the decays into massless vector bosons are all dominated by the 
Wess-Zumino contributions, which are far larger than those coming from the fermion loops. These 
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results will be used in the study of the relic densities of this particle which will be presented in the 
next sections. Here we compare the results of the decay rates for the new axions with those of the the 
PQ axion that we are going to compute from scratch. 

The interaction of the PQ axion with photons is given by 

where we denote with a the axion, which is bounded (from astrophysical and cosmological constraints) 
to be between 10^ GeV < fa "^ 10^^ GeV. The dots in the previous formula indicate terms that are 
irrelevant for the current analysis. From a general point of view, the coefficient c^^ depends upon the 
Peccei-Quinn charge assignment and also on the quark-mass ratios induced by its fermion interactions. 
The coupling to the fermions is given by 

TYl f — 

^f = %^V^/7'V'/, (30) 



where nif is the mass of the fermion, whose flavor is denoted by /, and the coupling gf = Qfj^ — Qf 



is given in terms of the chiral PQ charges {Q fi n) of each fermion (/). We denote with vpq the PQ 



breaking scale, which can be taken approximately around 10^^ GeV. We recall that in the PQ case 
the corresponding Wess-Zumino interaction is given by 

Ga 

4 



-^«77 ~ — A — " F^'^Ft^i^ — —Gay^a E ■ B, (31) 



where E and B are the electric and magnetic fields respectively, and the coupling Ga'yy is the sum of 
a model dependent term and of a second term which depends on the ratio of the quark masses 

G...^£|(E<?r(or)^-|i^). (32) 

where the quark-mass ratio is z = mu/md, while the Q'^'^'s are the e.m. couplings of the photons 
to the quarks. Since the coefficient Ga-yy is model dependent, we can have several possibilities. We 
compute below the decay rate into two photons in one specific case in which we assume 

47 = E^r(^r)' = o (33) 

/ 

and z = 0.56. This choice gives as a decay rate into two photons 

^.,, = |^^3^l.lxlo-^s-(^)^ (34) 

which is a function of the axion mass nia- More generally, we want to write the decay rate separating 
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Figure 2: Contributions to the x ~^ 77 decay. 



the contribution from the Wess-Zumino interactions from those which are obtained from the loop 
corrections. We obtain 



rpQ(a -^ -f-f) 



J2spin\-^PQ\ dki dk2 



2ma (27r)3A;0 {2TT)^k^ 

where the squared amphtude is given by 

/ ^ I-A^PqI = 2^ \Mpoint-like + Mloop\ 



(27rf6^^\k-ki-k2), 



(35) 



spm spin 



77 \ I ^ \ 4 , -"- 

3^ "^^+2 



F« 



?~^</"^'=«K-^) 



+ interf . , 



(36) 



where, in the second term, Nc{f) is the color factor, and the function r/ /(tj) is a function of the 
mass of the fermions circulating in the loop. We have introduced the function /(r), defined in any 
kinematic domain, whose real part is given by 



Re[/(T)] = I 
while its imaginary part is 

Im[/(r)] 



(arcsin 1/ y/rY 



i°«5' fe^ 



vr 



ifr > 1 
ifr < 1 



(37) 



M\m 



ifr > 1 
ifr < 1 



(38) 



where r = 4m^/m?,. In our case we take the branch r > 1. 

As we move to compute the decay of x ^-iid assume a free varying mass for this particle, the WZ 
interaction (Fig. [2^) is given by 



^wzix -^ 77) = 45'^7£[^, i^, ki,k2\. 



(39) 
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Figure 3: Total decay rate of the axi-Higgs for several mass values. Here, for the PQ axion, we have chosen 
fa = 101° GeV. 



In Fig. [2^ we have isolated the massless contribution to the decay rate coming from the WZ counterterm 
xF^F^ whose expression is 

(40) 



m" 



Combining also in this case the tree level decay with the 1-loop amplitude, we obtain for x ^ 11 
the amplitude 



M'^''{x^ll)=M^^z + M>}\ 
shown in Fig. [2j In this case the rates are derived from the expression 



(41) 



r(x -^ 77) 



327r 



8(5: 



11 > 



+ 



^-^ ATr'^nif 



^-l^Me'QJc^^f 



^^.^E^c(/).™e^g 



V fiv) ,2n2^xJ 



(42) 



and are shown in FigjSl In the equation above both the direct (~ (577)^) and the interference (~ g^-y) 
contributions are suppressed as inverse powers of the Stiickelberg mass. We show the results of this 
comparative study in Fig. [3l where in the left panel we present results for the decay rates of x — ^ 77 
for several values of the axion mass as a function of tan/3 = Vu/vd- The plots indicate a very mild 
dependence of the rates on this parameter, even for rather large variations. In the same plot the rates 
for the PQ case are shown as constant lines, just for comparison. Notice that we have chosen a rather 
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low Stiickelberg mass, with M = 1 TeV. The charge assignment of the anomalous model have been 
denoted as /(—1, 1,4), where we have used the convention 



These depend only upon the three free parameters qp , qf,Aq^. The parametric solution of the 
anomaly equations of the model /{qn ,qf,Aq^), for the particular choice qn = —l,qi = —1, 
reproduces the entire charge assignment of a special class of intersecting brane models (see 2j] and 



28l | and the discussion in jlSl ]) 

/(-I, -1, 4) = (-1, 0, 0, -1, 0, +2, -2). (44) 

In Fig. [3] (right panel) we show the decay rates as a function of the axion mass in both cases, having 
chosen a nominal mass range for this particle varying between 10~^ — 1 eV. One can immediately 
observe that the rates for the PQ case are smaller than those for the Stiickelberg by a factor of 
Iq'^o — 10^^, nevertheless the axi-Higgs x has a lifetime which is much bigger than the current age of 
the universe. 

Concerning the possibility to detect the axion through its two-photon decay channel, its tiny mass 
and the smaller value of its lifetime unfortunately do not allow to set significant constraints on its 
possible parameter space. The situation, in this case, if rather different from that of other dark 
matter candidates, such as, for instance, the gravitinos, which have been widely investigated recently 
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35l . 136( 1 . In fact, the allowed parameter space where the constraints derived from those previous 
studies apply, concern a region in the plane {T£)M,rn£iM) - with tdm being the lifetime of a generic 
dark matter particle and tudm its mass - which is bounded by the intervals 10^^ s < tdm < 10^^ s 
and 10~^ GeV < tudm < 10^ GeV. While the value of t^ for the axion can reasonably reach the lower 
edge of the scanned region in tdm, by an adjustment of its coupling gs and charge assignments of 
the anomalous C/(l), its mass is definitely too small to be excluded by these types of analysis. These 
studies are, obviously, very interesting for candidates of heavier mass, such as gravitinos. Similar 
considerations apply in the case of LHC studies, given the small production rates for a very light 
axion. For much heavier axions, instead, these types of studies have been performed quite recently 
[16t |. but the behaviour of this particle, in this case, is akin a light Higgs rather than a long-lived light 
pseudoscalar. 

5 Relic density at the electroweak and at the QCD phase transitions 

In this section we proceed with the derivation of the relic densities for x both at the electroweak and 
at the QCD phase transitions. 

At the electroweak scale, we will assume that the flat direction parameterized by the Stiickelberg 
axion b is lifted by electroweak instanton corrections. A similar phenomenon, but much more sizeable, 
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clearly will take place at the QCD phase transition. As we have discussed previously, at the electroweak 
scale, a mixing between the various phases of the non-perturbative potential allows to identify the 
linear combination x ^^ the physical axion. In general, this is misaligned with respect to the minimum 
of the potential generated at this transition, with a misalignment that, as we have pointed out, is 
parameterized by the value 9 = xl^x- 

The analysis of the relic density around the electroweak scale is then performed rather straight- 
forwardly, following a standard approach borrowed from the PQ case. For this goal, we define the 
abundance variable of x 

Tl 

Y^m) EE ^ (45) 

where Tj is the oscillation temperature, which is close to the electroweak scale. The universe must be 
(at least) as old as the required period of oscillation in order for the axion field to start oscillating and 
to appear as dark matter, otherwise 6 is misaligned but frozen. This is the content of the condition 

m^{Ti) = 3H{Ti), (46) 

between the mass of axion at the oscillation temperature Tj {my.{Ti)), and the Hubble parameter at 
the same temperature H{Ti). The condition for oscillation Eq. ()46|) allows to express the axion mass 
at T = Tj in terms of the effective massless degrees of freedom evaluated at the same temperature, 
that is 

-x^) = \lr'9.,T.§^. (47) 

Expressed in terms of the initial angle of misalignment 9i , Eq. [35] becomes 

YJTi) = ^ ' , (48) 

where g-t^T = 110.75 is the number of massless degrees of freedom of the model at the electroweak 
scale. Using the conservation of the abundance YaO = Ya{Ti), the expression of the contribution to 
the relic density is given by 

T, Pc 

The values of the critical energy density (pc) and the entropy density today are estimated as 

Pc = 5.2 • W-^GeV/cm^ sq = 2970 cm-^ (50) 

with ~ 1. Given these values, the relic density as a function of tan/3 is given in Fig. [H We have 
varied the oscillation mass and plotted the relic densities as a function of tan (3. The variation of Vu has 
been constrained to give the values of the masses of the electroweak gauge bosons, via an appropriate 
choice of tan/3. 
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Figure 4: Relic density of the axi-Higgs as a function of tan/3 

For instance, if we assume a temperature of oscillation of Tj = 100 GeV, an upper bound for the 
axi-Higgs mass, which allows the oscillations to take place, is ?7i^(Tj) ~ 10~^eV, with g^^x ^ 100. 

In order to specify a^. we have assumed a value of 1 TeV for the Stiickelberg mass M, with 
gs ~ 1, and we have taken (g„, q^) of order unity, obtaining o"^ ~ 10^ GeV. As we lower the oscillation 
temperature (and hence the mass), the corresponding curves for Q^ are down-shifted. 

The values of these relic densities at current time are basically vanishing and these small results 
are to be attributed to the value of a^, which is bound to vary around the electroweak scale. 

Just to compare with the PQ case, there a^ is replaced by the large scale fa at the QCD phase 
transition, and this is the reason of such a strong suppression for i}^ (or of an enhancement, in the PQ 
case). Instanton effects at the electroweak scale are expected, in our case, to provide a mass of the type 
m^ ~ Ag^/v^, with Ag^ ~ Exjp{—2TT/aw{v))v'^ - aw{v) being the weak charge at the scale v - which 



is indeed a rather small value since Exp(— 27r/au,(w)) 



~ e 
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. For this reason x remains essentially 



a physical but frozen degree of freedom which may undergo a significant (second) misalignment only 
at the QCD phase transition. If not for the presence of a coupling of the axion to the gluons, via the 
color/ U{1)b mixed anomaly, x could be classified as a quintessence axion, contributing to the dark 
energy content. 

5.1 The QCD phase transition 

We have seen that the electroweak phase transition has redefined the phase of x via the mixing with the 
CP-odd Higgs sector, but below the electroweak scale the field remains essentially a pseudo Nambu- 
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Goldstone mode which undergoes the second niisahgnment induced by the QCD phase transition, 
quite similarly to an ordinary PQ axion. Neglecting the small mass of m^ induced at the electroweak 
scale, the new mass induced at the QCD scale is controlled by the ratio m^ ~ A^^^^u/M^, where now 
the angle of misalignment is essentially related to the Stiickelberg mass M, via M'^ jv (which replaces 
cj;^), and is now given by 0' = x^/-^^- 

In order to further clarify this point it is convenient to follow the analogy with the PQ case and 
observe that M^/v replaces /„ in characterizing the coupling of the physical axion x to the gluons. 
At the same time M, just like /a, can be interpreted as a symmetry breaking scale, given the presence 
of a derivative coupling {Mdh ■ B) of the anomalous gauge boson B to the Stiickelberg field b in the 
Stiickeberg mass term. Thus it can be naturally interpreted, in this phase, as originating from a vev 
of an extra scalar singlet to which B couples in the UV. These two elements clearly indicate that the 
new misalignment is basically given by 9' . 

Given the similarity between this situation and the PQ case, then we can follow standard arguments 



to estimate the mass of x after the QCD transition. Thus, we just recall that for PQ axions [3|, [37 1 
the zero temperature mass is given by 

ma = N-— — = 6.2- /xeV (51) 

^ + Z fa Ja,12 

where z ~ 0.56 is the ratio of the up and down quark masses, /jr and m^ are the pion decay constant 
and mass, fa is the axion decay constant, /a, 12 is the same constant expressed in units of 10^^ GeV and 
N is the U{l)pQ color anomaly index. The dependence of the Peccei-Quinn axion on the temperature 
can be expressed as 

[ma r < A, 

where d is a model dependent numerical factor and Aqcd ~ 0.2 GeV is the scale of the QCD phase 
transition. We have also set d = 0.018 [38]. 

We can borrow this formula to determine the mass of x at zero temperature {m^{T = 0) ~ m^) 
and extend it to finite T (m^{T) ~ ma{T)), using the same expression ([52]) valid in the case of the PQ 
axion. Coming to the value of the abundances, with the replacement of a — )■ M^/f, Eq. HHl assumes 
the form 

YJTi) = — ^ / ^ ^ ^ -^, 53 

being the exact analogue of the PQ expression for the abundances of the invisible axion, with fa — >• 
M"^ /v. Concerning g^^^x, the effective massless degrees of freedom at T ~ 1 GeV are those of the gluons, 
the photon, 2 charged leptons, 3 neutrinos and 3 quark flavors; thus we have g^^^iceV = 61.75. At the 
QCD phase transition, that is at T ~ 200 MeV, the effective massless degrees of freedom are given by 
the photon, 1 charged leptons and 3 neutrinos, giving g^^\ = 10.75. The oscillation temperature can 
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be obtained from Eq. UT] 



Tf = V 4^^4:^"'x(0)Mp hhtiCD, (54) 



from which we get Tj ~ 0.6 GeV and a typical oscillation mass which is given by m^(Ti) ~ 1.4 neV. 
Values of m^ larger than this typical value will allow oscillations of the field x ^^d the appearance 
of relic densities whose size is essentially controlled by the value of M, the Stiickelberg mass, via the 
scale M'^/v. 

Given the analogy between M'^/v and fa and the dependence of the axion field amplitudes on 
these two scales, it is natural to expect that only for large values of M one should expect a significant 
contribution to the relic density of these new axions. 

We show in Fig. [5] results of a numerical study of ^mish? as a function of M, expressed in units of 



10^ GeV. We show as a darkened area the bound coming from WMAP data 39|], given as the average 
value plus an error band, while the monotonic curve denotes the values of Qmish"^ as a function of M. 
It is clear that the relic density of x can contribute significantly to the dark matter content only if 
the Stiickelberg scale is rather large (~ lO'^ GeV) and negligible otherwise. A final comment concerns 
the role of the isocurvature perturbations, which are generated by infiation, in these types of models, 
since in the case of the PQ axion they provide significant constraints on the possible values of fa- The 
fact that the b field does not correspond to a physical degree of freedom during inflation allows to 
bypass completely these constraints. They do not apply to these types of axions and this represents 
a very interesting feature and a significant variant of these models respect to the PQ case. 

6 Conclusions 

We have discussed the most salient cosmological features of models containing gauged axions, obtained 
from the gauging of an anomalous symmetry. The gauging allows to define a consistent theory for 
axion-like particles, which generalize many of the properties of PQ axions. They have appeared for 
the first time in the study of intersecting branes, but their features are quite generic. They are 
constructed as effective theories containing minimal gauge interactions which restore gauge invariance 
of the effective action in the presence of an anomalous f/(l) symmetry, and no further requirements. 
Differently from the PQ case, here there is no concept of an original PQ symmetry, broken at a 
very large scale, with the axion taking the role of a Goldstone mode that acquires a mass at the QCD 
phase transition. Rather, the physical axion emerges directly at the electroweak phase transition, when 
Higgs-axion mixing occurs. Being charged under SU{3) and SU{2), we have a sequential misalignment 
of this field, and we have quantified its relic density as a function of the Stiickelberg mass. We have 
shown that only very large values of the Stiickelberg mass cause a significant contribution of this type 
of axions to the current dark matter content of the Universe, which otherwise remains negligible. The 
absence of an original PQ-like potential has some implications at cosmological level, such as the absence 
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Figure 5: Relic density of the axi-Higgs as a function of M. The grey bar represents the measured value of 
^DAih^ =0.1123 ±0.0035 

of isocurvature perturbations, since the Stiickelberg is not a physical mode before the electroweak phase 
transition, in particular at the time of inflation. This feature is due to the presence of a local gauge 
symmetry, realized in the Stiickelberg form, which allows to absorb b into the longitudinal component 
of the anomalous gauge boson. 

Our analysis represents, more generally, a description of the fate of the Stiickelberg field in cos- 
mology, from the defining Stiickelberg phase of the theory at a large scale (defined by the value of 
the Stiickelberg mass), down to the electroweak and QCD phase transitions, when this field devel- 
ops a physical component. Our analysis could be extended in several directions, for instance with 
the inclusion of the modifications induced on the computation of the relics due to the presence of 
non-homogeneities in x beyond the QCD horizon, a feature which is also present in PQ models when 
the PQ scale lays below the scale of inflation. However, even at this level of refinement, the only 
significant scale in the determination of the relic densities remains the value of the Stiickelberg mass. 
Small values of this mass parameter in the TeV range leave the contribution of these particles to the 
relic densities of dark matter negligible, and sizeable for M around an intermediate scale of 10^ GeV. 
In this case all the constraints coming from the neutral current sector are satisfied, being the extra 
Z' of the theory completely decoupled from the low energy spectrum of the Standard Model. The 
appearance of this intermediate scale is a novel feature of this type of axions which could be used to 
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set limits on their parameter space. 
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A Appendix. The model, definitions and conventions 

We summarize in this section some results concerning the model with a single anomalous U{1) discussed 
in the main sections. 

The effective action has the structure given by 

S = Sq + Svuk + San + Sw Z + Scs (55) 

where Sq is the classical action. It contains the usual gauge degrees of freedom of the Standard Model 
plus the extra anomalous gauge boson B which is already massive, before electroweak symmetry 



breaking, via a Stiickelberg mass term, reviewed in Sec. [2l Its complete expression is given in 171 ]. 
Here we briefly describe the structure of the anomalous contributions and of the induced counterterms 
for the restoration of gauge invariance in the 1-loop effective action. 

In Eq. ()55p the anomalous contributions coming from the 1-loop triangle diagrams involving abelian 
and non-abelian gauge interactions are summarized by the expression 

San = ^^{TbwwBWW) + ^^{TbggBGG) + j^{TbbbBBB) 

+^^{TbyyBYY) + ^^{TybbYBB), (56) 

where the symbols () denote integration. For instance, the contributions in configuration space are 
given explicitly by 

{TbwwBWW) = JdxdydzTX§{z,x,y)B^{z)Wt{x)W;'{y) (57) 

and so on, where Tbww denotes the anomalous triangle diagram with one B field and two VF's external 
gauge lines. The gluons are denoted by G. The Wess-Zumino (WZ) counterterms are given by 

SwZ = CBB{bFBAFB)+CYY{bFY AFY)+GYB{bFY AFb) 

+F{bTr[F^ A F^]) + D{bTr[F^ A F^]), (58) 

while the gauge dependent Chern-Simons (CS) abelian and non abelian counterterms [4fl] needed to 
cancel the mixed anomalies involving a B line with any other gauge interaction of the SM take the 
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form 



Scs = +di{BYAFY)+d2{YBAFB) 



The non-abelian CS forms given by 

(^SU{3) ^ 1 

• The structure of g^^ 



K (Cp + ^92 e'^'WlW^^ + cyclic 



T^G I „ fahc/-~ih /^c 



Gl[F^^,p + -g^nGlGl]+ cyclic 



(60) 
(61) 



The coefficients in front of the WZ counterterms are determined by requiring gauge invariance of the 
effective action. We outhne the case of g^-y and its relation to the fundamental parameters/scales of 
the theory. Among these are the Stiickelberg mass M, the hypercharge and weak couplings gy and 
52 and the charges of the fermion running inside the anomaly loops. These fix the coefficient of the 
anomalies Cbyy and Cbww (for the U{1)b U{1)\ and SU{2)'^ U{1)b anomalies) and the rotation 
matrices of the neutral gauge bosons O and of the CP-odd sector O^, defined in Eq. (J16p . This is 
defined as in Eq. ()20|] in terms of the counterterms 

F = ^igl^CBww, (62) 

with 

Cbww = ~^z_^'ifL-, (63) 

/ 

with Qn = —-^ being the AW anomaly, and 

Cyy = ^-^igf-^CBYY. (64) 

which is defined by the charges 

Cbyy = \ E {<llMnf " ^U^\f\ ■ ^ 

f 

Explicit expressions for Cbyy and Cbww are given in Eq. ([9]). 
• Fermion interactions 

The covariant derivatives are defined as 

D^ = d^ + igsT'^Gl + i52T'^T^; + '-gYQ^Y^ + '-gBQ^'B^, (66) 
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with T^ and r^ given by 



rpa 



x° 



(67) 



2 2 ' 

where A" and a" are the Gell-Mann and Pauh matrices. This choice of the covariant derivative defines 
the gauge variations of the fields; in particular, under the abelian group transformations we have 



B'^ = B^ + d^e b' = b + Me ^' = e-'^ss?.^ ' 

We write the lepton doublet as 



(68) 



L, 






The interaction Lagrangian for the leptons becomes 



C 



lep 

int 



"1 



-B , 



-92t''W; + -qyY^ - -gBqfB^ 






+ 



ena^ 



,B 



j5r^M - ^gsqei^Bf, 



e-Ri- 



As usual we define the left-handed and right-handed currents 



J, 



r)\'^fJ. '^IJl)i 



•^^ ~ 9 \'^^^ + "^ fi) 






r - T" 



J' 



Writing the quark doublet as 



Q^ 






we obtain the interaction Lagrangian 



C 



quarks 
int 



ULi dLi]l^ 

+ una'' 



-gsT'^Gl - g2T-W^ - -gyY^ - -gsq^B^ 



-|2-^^ -p ryJiJru^-'lJ- 



ULi 

.dii 



+ 



-gsT-Gl - g2T-W^ + ^gyY^ - \gBqlB, 



URi 
dRi- 



(69) 



(70) 



(71) 



(72) 



(73) 



We work with a 2-Higgs doublet model, and therefore we parametrize the Higgs fields in terms of 8 
real degrees of freedom as 



Hu 



Bu 



H, 



Hi 



(74) 
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where H^, H^ and i7„, if^ are complex fields. Specifically 

+ _ Reff + + ilmH+ _ Reff' + ilmff- 

~ V2 ' '^ ~ 71 ' ^u - ^u 1 ^d - ^d ■ I'^i 

Expanding around the vacuum we get for the uncharged components 

,_ Reffg + ^Imffg ,_ E^Hl+JJmHl 

^u-Vu-\ y= , H^-Vd^ -7= . (Tbj 

The Weinberg angle is defined via cos^vi^ = 92/ 9-,sin9w = gv/g-, with g"^ = Qy + si- ^^ ^^^^ define 
cos/3 = frf/w, sin/3 = Vu/v with w^ = w ^ + w^. 

A.l The Yukawa couplings and the axi-Higgs 

The couplings of the two Higgs and of the axi-Higgs to the fermion sector are entirely described by 
the Yukawa Lagrangian. The Yukawa couplings of the model are given by 

^urut. ^ _^d Qjj^^^ _ pd J^^tg _ r« QL{ia2Hl)uR - F" UR{ia2H*jQL 

-T' LHden - V ErHIl 
= -T'^ dH^^PRd - V^ dHfPLd - F" uHI*Pru - F" uHIPlu 

-T' eH^PRC - r eH^/PLe, (77) 

where the Yukawa coupling constants F , F" and F*^ run over the three generations, i.e. u = {u,c,t}, 
d = {d,s,b} and e = {e, n, r}. Rotating the CP-odd and CP-even neutral sectors into the mass 
eigenstates and expanding around the vacuum we obtain 

ReH^ + UmHl 

Hu=vu + ^ 

, {h^ sin g - g" cos a) + i {O^.Gl + 0^01 + Ofg) .„«. 

-Vu+ ^ U«j 

H'd =Vd + ^ ^ 

^^^ ^ (/i° cos « + ff sin a) + i {O^^Gl + O^^Gl + O^^x) ,^^. 

\/2 



so that in the unitary gauge we obtain 



ijO =vu + —1= [(/i° sin a- H^ cos a) + i Ofix] 
v2 

H^d =Vd + -^ [(/i° cos Q + F° sin q) + i0^2x] , (80) 

v2 
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where the vevs of the two neutral Higgs bosons Vu = v sin /3 and Vd = v cos /? satisfy 

tan/3 = ^, v = Jvl + vl (81) 

Vd ^ 

The ferniion masses are given by 

m„ = VuV^, ruu = VyT" , 

ma = VdV'^, me = VdT', (82) 

where the generation index has been suppressed for brevity. The fermion masses, defined in terms of the 
two expectation values v^, Vd of the model, show an enhancement of the down- type Yukawa couplings 
for large values of tan /3 while at the same time the up-type Yukawa couplings get a suppression. The 
couplings of the h^ boson to fermions are given by 

cosa,o\ ™„-, ,. /sma,o\ ™, „ /cosa.Q 



>CYuk(/i") = -r-^ dLdR ( ^T^h'] - r" ULUR {-^h') - r" ^lsr i-y^h" ] + c.c. (83) 



The couplings of the H boson to the fermions are 



rO^ _ rdj , /sma o\ ^u„-.„. / cosa o\ r-e ^ ^ /smaj^o 



CYni.{H") = -TUlcIr ( -^H") - ^""ULUR ( -— =-FU 1 - T'eLCR {-^H'> ) + c.c. (84) 



The physical gauge fields can be obtained from the rotation matrix O 



A 




which can be approximated at the first order as 




(85) 




where 



9y_ 92 

0^:.|f +0(ef) -^ + 0(6?) §61 I (86) 



_ XB 

XB = {q^vl + q^vj) . (87) 



More details can be found in 



Q. 
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